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Abstract —This paper puts forward the vision of fully decou¬ 
pling market operations from security considerations through 
controllable power flows. In “A Fully Controllable Power Sys¬ 
tem”, power system security is no longer dependent on the 
location of the power injection points. In the ideal case, this leads 
to the elimination of redispatching costs, which amount to several 
million dollars per year in large systems. This paper determines 
the upper and lower bounds for the number of controllable 
lines and number of controllers to achieve this decoupling in 
any system. It further introduces the notion of the controllability 
vector CV, which expresses the effect of any controller on the AC 
line flows. Based on two alternative deflnitions for controllability, 
two controller placement algorithms to maximize controllability 
are presented and their results are compared. 

Index Terms —High-Voltage DC lines, FACTS, controllability, 
optimal power flow, optimal placement 

1. Introduction 

“A Fully Controllable Power System” aims to arrive at 
a system where the market operations are fully decoupled 
from the security considerations. In such a system, the “Cost 
of Security” will be zero, as the system controllability will 
guarantee the security of the system without need for any 
costly preventive control actions, such as the generation re¬ 
dispatching. 

The first objective in power system operation is to satisfy the 
power demand. As long as a Standard OPF results in a feasible 
solution, the system has both the necessary generation and 
transmission capacity to supply all loads. Normally, however, 
most power systems are dispatched in a way to be N-1 secure. 
For that purpose, a Security-Constrained OPF (SCOPF) is 
executed. Due to the uncontrollable nature of the AC line 
fiows, the points where the power is injected play a crucial 
role on the resulting power fiows and on possible overloadings. 
As long as the Security-Constrained OPF results in a feasible 
solution, the system has additional generation capacity, and, 
most importantly, sufficient transmission capacity, so that even 
if an element is lost, no overloadings will occur. By definition, 
the SCOPF will result in an equal or more expensive dispatch 
than the Standard OPF. In case we were able to control all line 
fiows, however, the injection points would play a substantially 
less significant role. 

The “Cost of Security” reflects the additional costs incurred 
to the system, so that N-1 security is ensured Q. It is 
calculated as the difference between the system costs resulting 
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from a Security-Constrained OPF minus the system costs 
resulting from a standard OPF, as shown in ([T]). 

CoS = CsCOPF — CoPF (1) 

Several approaches exist in the literature, which aim to take 
advantage of power system controllability to provide both pre¬ 
ventive and corrective control and minimize the redispatching 
costs. Ref. n presented the first formulation incorporating 
corrective control actions from generation rescheduling and 
switching actions. In 13 the benefits of transmission switching 
for optimizing the power system economic operation are 
demonstrated. In n optimization algorithms for corrective 
control with FACTS are presented, while, to our knowledge, 
m presented the first SCOPF formulation with corrective 
control of HVDC lines. 

The goal of the concept “A fully Controllable Power Sys¬ 
tem” is to decouple the dependence of the line fiows on the 
power injection points. Market operations, e.g. the solution of a 
standard OPF, should determine the most economical dispatch. 
Through fully controllable power flows, system operators will 
guarantee system security by appropriately rerouting the power 
before and/or after a contingency occurs. Thus, as long as the 
transmission capacity is sufficient, there will be no need for 
(preventive) generator redispatching actions and the dispatch 
determined through the market is accepted as is. By incuring 
no additional redispatching costs this results to a zero “Cost 
of Security” in the ideal case. 

In the rest of this paper we will deal with two main 
problems. First, the minimum number of controllers that are 
necessary in order to make the system fully controllable will be 
determined. Second, placement algorithms in order to achieve 
maximum controllability in the system for a given number of 
controllers will be introduced. 

The following methods focus mostly on HVDC lines. How¬ 
ever, with small modifications they can also be applied to 
Thyristor-Controlled Series Capacitors (TCSCs) and Phase- 
Shifting Transformers (PSTs), as we will discuss later. In 
the text, we will use the terms “HVDC” and “controllers” 
interchangeably. 

H. Deeinition of Controllability 

In this work, we aim at maximizing the controllability of 
the system for a given number of controllable elements. We 
define controllability as the magnitude of the change in line 
fiows caused by a marginal change in the controller setpoint. In 
order to maximize controllability in the system we are seeking 
controller placements that will have the largest impact on line 
fiows and require the smallest control action. 
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In mathematical terms, for a single controller j we define 
controllability maximum absolute 

change on each AC line flow for control actions 

I APc(j) I < P where i G £ = {1,... ,72^}: 

riL 

P(i) = |APc(i)|<P- (2) 

i=l 

Based on (|2]), we are trying to maximize order 

to maximize controllability. If the system is linear (e.g. DC 
approximation), ^c(j) linear in p. 

III. The Controllability Vector CV 

For the derivations in the following sections we will follow 
the DC approach (i.e., linearized equations for power flows) 
and use the Power Transfer Distribution Factors (PTDFs). 
Nevertheless, we can derive very similar relationships using 
the full AC approach, as shown in O. For that we need to 
focus on the bus injections currents and the line currents, and 
make use of the Kirchhoff current law. 

The PTDFs are linear sensitivities, which express the effect 
of bus injections on the line flows, as shown in dS]). PTDF 
is an ul x ub, where is the number of lines and ub the 
number of buses H, Q. 

Pl = PTDF Pb (3) 

The equations for the HVDC links follow the formulations 
introduced in (8) and 0. Each HVDC link is modeled as two 
virtual voltage sources located at the two nodes where the 
HVDC line is connected. For each HVDC link, we assume 
one additional variable for the HVDC power flow. Let Pdc ^ 
]^noc represent the vector of power flows on the HVDC lines. 
For an HVDC link j connected between nodes m and n, the 
balance between the active power injected or withdrawn from 
the line is maintained by assuming that Pinj,m = —Pinj,n = 
PdcU) (consistent with the DC power flow approach, HVDC 
line losses are here neglected). 

For an arbitrary change APdc in the power flow of the 
HVDC line, the following equation will hold: 

APl = (PTDF„, - PTDFn) APdc. (4) 

where PTDFm, PTDFn are x 1 vectors corresponding 
to the columns m and n of the PTDF matrix. 

Based on 0, we define the controllability vector CV to 
represent the effect of a marginal change of the HVDC line 
flow on each AC line of the system, as follows: 

CV„,n = PTDF„, - PTDFn (5) 

Each element of the controllability vector represents the 
change in the power flow of a specific AC line resulting from 
the change of the HVDC line setpoint. The higher the absolute 
value of this element is, the larger the influence that the DC 
line has on this specific AC line. In Section [2 we also derive 
the controllability vector for TCSCs. 

There exists one controllability vector for each pair of nodes 
in the system. As a result, each HVDC placement between two 
arbitrary nodes corresponds to a unique controllability vector. 
The controllability vector will become central in our analysis 


and the design of the placement algorithm in the following 
sections. 

In the following sections we will derive the upper bound 
for the number of controllable lines in a system. We will 
distinguish between the placement of controllers in series and 
in parallel. Series controllers could be TCSC devices, HVDC 
back-to-back converters, or replacing an existing AC line with 
an HVDC line. Controllers in parallel are new HVDC lines 
that either connect two new pair of nodes or are placed in 
parallel to an existing line. 

IV. Maximum number of controllable lines for 
FULL controllability: Placement of Series 
Controllers 

Power grids have a planar graph topology. It has been shown 
that any planar graph can be decomposed to a set of series, 
parallel, and Wheatstone graphs ifTOl . Figure [H shows the three 
different graph types. Showing that the following derivations 
apply to each of these graph types, it follows that it will apply 
to any power grid topology which is a synthesis of such graphs. 



Fig. 1. Graph decomposition. Every graph can be decomposed to series (left), 
parallel (middle), and wheatstone (right) graphs. 

We assume a power system with n b nodes and ul branches 
connecting the nodes. For every node k it will hold: 

Pk = Pki + Pkj + Pkm + • • • (6) 

where j, /c, m,... G {1, ... ,72^}. Pki. Pkj. • • • represent the 
line flows on lines connected at node k, while Pk represents 
the net power injections (resp. withdrawals) at the same node. 
If no shunt element (e.g. generator, load, shunt capacitance, 
etc.) is connected to a node k, then Pk = 0. 

We now assume that the system is fully controllable. This 
means that given a certain profile of power injections at all 
nodes, we can determine independently the power flows on 
each line. Such a system could be thought of as a network of 
HVDC point-to-point links. 

The algebraic system of (O V/c G 72^ is of the form Ax = 
B: 

Pb = APl (7) 

Aisa72BX72L matrix whose elements are ones or zeros. Pb 
is a 72 b X 1 vector equal to the bus power injections: Pb = 
[Pl... PubV- Pl corresponds to a 72^ x 1 vector containing 
all the line flows, i.e. Pl = [Pij], with iC G {1,... ,72b}. 

We assume the net injections/withdrawals are given, so 
Pb is a vector of constants. The question is how many 
elements of the Pl vector should be known (or “exogenously 
constrained”) in order to be able to fully define the power 
flows, i.e. the rest elements of the Pl vector. 
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For a system of equations to have a unique solution, matrix 
A must be square and non-singular. In a meshed power system, 
however, there are usually more branches than nodes, i.e. > 
nB . This means that our problem is underdetermined. 

Assuming no losses in the lines, for all nodal power 
injections Pk must hold ^j^Pk = 0. This implies that the 
rows of A in © are linearly dependent. In total, we have 
ub — I linearly independent equations. We assume here that 
all parallel lines between two nodes are represented by one 
equivalent line. Thus, matrix A is of full column rank. 

If riL < riB At must hold ul = ub — f, since each branch 
connects two nodes. Having also ub — f linearly independent 
rows, matrix A is square and non-singular, and our problem 
yields one unique solution. Consequently, for a given set of 
bus injections, ifn^ = ub — I (e.g. see series graph in Fig.[T]) 
none of the line flows can be controlled to an arbitrary setpoint. 

If riL = riB and we have ub — f linearly independent 
equations, we need to deflne one more equation in order to 
have a unique solution of the system. As a result, if = ub 
we need one controller to set one line flow to a speciflc value. 
According to this value the rest of the line flows can be 
calculated. 

If > Ub, then we can freely select +1 variables. 

As a result, even in a fully controllable power system we can 
freely determine the flow only on ul — tib + f lines. The rest 
of the line flows will be determined through the bus power 
injections (see also the parallel and Wheatstone graphs in 
Fig.©. This implies that during a controller placement process 
- converting an AC non-controllable system to a system with 
controllable power flows - we would not need to place more 
than riL — riB P f series controllers to control the maximum 
number of lines for a given set of bus power injections. 

This should hold, as long as the controllers are assumed 
to have inflnite (i.e. sufficiently large) boundaries and there is 
sufficient transmission capacity. If their operation is limited, 
either we might need more controllers, or we might never 
achieve full controllability of the power flows. This leads us 
to the following two conclusions, as shown in ([5]). First, by 
placing controllers in series with existing AC lines, we cannot 
control more than ul — tib P f branches. Second, in order to 
achieve full controllability we need at least nL — riBpf series 
controllers: 

= nL-nB + l = inf A''conttollers- (8) 


as shown in I©. 


PL(fc + l) = I-PL(fc) + CV-PDc(fc) (9) 


Then the controllability matrix C is equal to: 

C = [CV I ■ CV ... F^-^CV] (10) 


From ([Toll , it is rank{C) = rank{C'V). However, CV is only 
a linear transformation of PTDF, i.e. CV = T • PTDF, 
where T has dimensions x noc and is of full rank, since 
it is a set of linearly independent vectors containing 0,1, — 1. 
Assuming that there are controllers installed to control flows 
between every possible pair of nodes, then udc = TiB{nB — 
l)/2 > riB and rank{CY) = rank(PT'DF) = — 1. 

This means that for a given set of power injections, we can 
control up to — 1 lines by installing parallel power flow 
controllers, see ([TT]) . For example, in the series graph of Fig. [T] 
we can install up to four parallel controllers and we will be 
able to control different degrees of freedom of the system. 


TVT-parallel -• • n Tvrparallel 

S«P^ctrl.lines = «£?-!= inf 


controllers' 


( 11 ) 


As we will see in the following sections, each controller 
(either DC lines or TCSC) can only control a single line. As 
a result, even in the ideal case where all line flows were fully 
controllable and the controllers had inflnite bounds, at any 
given point no more than ub — f controllers could be active. 


VI. Adding a controller to the system 

In this section we now assume that we have a non- 
controllable AC network. We further assume that the line 
capacities are sufficiently high so that no congestion is induced 
by different bus power injections. We will focus on two 
main power flow control elements: the HVDC line and the 
Thyristor-Controlled Series Capacitor (TCSC). We will show 
that in the general case each of these elements can determine 
the flow on almost any single AC line, but only one at a time. 


A. HVDC line 

From 0 and (O, it follows that: 


V. Maximum number of controllable lines for 
FULL controllability: Placement of Controllers 
IN Parallel 

According to control systems theory, full controllability 
describes the ability of an external input to move the internal 
state of a system from any initial state to any other flnal state 
in a flnite time interval. We can determine how many states 
of the system are controllable from the controllability matrix. 
If the controllability matrix has full row rank, then the system 
is controllable. We deflne the discrete-time invariant systeirQ, 

^The results would be similar for a continuous LTI system. Then x = 
Ox(f) + Bu{t). The controllability matrix is C = [B] = [CV], and 
rank(C) = ub — t 


PB = CVmn^P^A ( 12 ) 

where m, n are the nodes connected through the HVDC line. 
With P^^ being a scalar, there is only one degree of freedom, 
and therefore each HVDC line can control only a single AC 
line at a time. Through (fT^ . we observe that by adding an 
HVDC line, we can influence all AC power flows in the system 
as long as all elements of the vector Cmn are non-zero. This 
is most often the case. In a meshed transmission system it is 
relatively rare that power injection and withdrawal between 
a pair of nodes m,n influences in exactly the same way a 
line flow k (so that PTDFm,k, PTDFm,k, of CVmn,k equal 
zero). As a result, by adding one HVDC line we can control 
the flow of almost any single AC line. 


4 


B. Thyristor-Controlled Series Capacitor (TCSC) 


A TCSC is connected in series to a transmission line and 
is able to control the reactance of this line. As a result, it can 
indirectly control the line flow either on this line, or influence 
the flow of any other line in the system. Through the following 
derivations we show that one TCSC device can influence 
almost all line flows in the system (to different degrees) and 
can arbitrarily determine the line flow on exactly one line. 

In a DC power flow context it holds: 


Pl = PTDF Pb Pl = BlB-^Pb (13) 

Bl is the line susceptance matrix. B^^ is the inverse of the 
bus susceptance matrix BbE|. 

We wish to examine the change in the line flows Pl, by 
changing the line reactance Xij. It is: 


^Pl 

dxij 


^Pl 

dxij 

^Pl 

dxij 


dxu 


+R R-i 

= 0 


+ BbB5> ^B5‘Pb 


dXi 


+ PTDF^^^Br^Pb 


dx. 


dxi 


(14) 


(15) 


(16) 


From ([T^ . the term PTDF^^^Bg^ will always result 
in a non-sparse matrix (in the general case where no radial 
connections exist all values will be non-zero, except for the 
column corresponding to the slack bus). This leads to the 
conclusion that in a meshed system a TCSC can influence 
almost all AC lines in the system. Given that one TCSC has 
a single degree of freedom, assuming that it has sufficiently 
large boundaries, it can determine arbitrarily the flow on any 
single line in the system. 



Fig. 2. 10-bus network used for the simulations. 


(CoS) again. We do that for each possible line placement, 
before we move on to the installation of two TCSC devices 
at the same time. All possible combinations of line pairs are 
examined. 

Figure [3] presents the Cost of Security in percentage, after 
placing one TCSC, two TCSCs, etc. up to 14 TCSC devices. 
The values shown are for the TCSC placement that achieved 
the highest CoS reduction from all possible combinations at 
each placement step. With the installation of TCSC devices 



VII. Example for Placement of Series 
Controllers - 10-bus network 

Consider the power system model of Fig. [51 It consists of 
10 nodes and 14 AC lines. The system data can be found 
in the Appendix of ifTTIl . In this example we will examine the 
placement of TCSC devices. We will show that, given a single 
generation and demand snapshot, there is a maximum number 
of controllable devices we can install in a power system, above 
which no additional power flow controllability can be gained. 
This upper bound is equal toni,—rfB + l. For this speciflc 
case study, with = 10 and = 14, this amounts exactly 
to 5 TCSC devices. 

We start by having no TCSC device installed in our system. 
We run a Standard AC-OPF and a Security-Constrained OPF 
and compute the Cost of Security - as deflned in ([T]) - in 
percentage values. For the SC-OPF algorithm, we use the 
formulation described in HD- Subsequently, we add a TCSC 
device along an existing line and compute the Cost of Security 

^Because the bus susceptance matrix is singular, before inverting it we 
eliminate the row and the column corresponding to the slack bus. After the 
inversion we insert zero vectors at this row and column which correspond to 
the slack bus. To denote this manipulation, we express the inverse of Bb 
with a tilde, i.e. 


Fig. 3. Cost of Security (in % relative to AC-OPF Costs) for different number 
of TCSC devices installed 

both the AC-OPF and the SC-OPF costs decrease. At the same 
time, their difference, both in absolute and percentage terms, 
also decreases. As shown in Fig. (3 after the placement of the 
first TCSC device we observe a significant reduction in the 
CoS. The controllability introduced through this TCSC offers 
additional flexibility in the routing of the power flows which 
reduces the amount of redispatching necessary to maintain 
the system N-1 secure. It is interesting to note that after the 
placement of five TCSC devices, the CoS stabilizes and no 
further significant reduction on the CoS can be observed, e.g. 
by adding six or more TCSCs. 

According to our considerations, the minimum number of 
TCSC devices for full controllability is 5 controllers (jil = 14, 
riB = 10, see ([8])). This is confirmed from the simulation 
results shown in Fig. [3 The placement of further controllers 
does not add anything to the controllability of this system, and 
it has no effect for the cost reduction of the speciflc snapshot. 

Two additional comments about this case study should 
be made here. First, as we can observe, the CoS is not 
eliminated but amounts to about 0.72% of the AC-OPF costs 
after the placement of five or more TCSCs. Second, the CoS 
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“stabilizes” at about 0.72% but small oscillations around this 
value can be observed with the placement of additional TCSC 
devices (i.e. above five). There are two reasons for these two 
phenomena. First, in order to completely eliminate the Cost 
of Security not only controllability plays a role but also the 
transmission capacity. In case of contingency, enough trans¬ 
mission capacity should be available, so that the power could 
be rerouted through alternative paths. This means that in case 
of higher loading in the system the CoS could “stabilize” at 
values higher than 0.72%, while for lower system loading CoS 
can equal zero. The second reason has to do with the TCSC 
limits. For this case study, we have selected realistic limits 
for the control capabilities TCSC devices: the TCSC can vary 
the line reactance up to 60%. With sufficiently large limits, 
the CoS could be completely “stabilized”, and eliminated if 
there is sufficient transmission capacity. An additional reason 
for these small “oscillations” is the fact that both the AC-OPF 
and the SC-OPF are non-linear and non-convex problems. As 
a result, there might be small differences due to numerical 
reasons, as the solver might not be able to find the absolute 
optimal solution. 

VIII. HVDC Placement eor Maximum 
Controllability: Placement of Parallel 
Controllers 

In the rest of this paper, we propose algorithms for the 
placement of HVDC lines to maximize the controllability in 
the system, based on two different definitions of controllability. 
The goal is to achieve the maximum controllability in the 
system with the minimum number of controller placements. 

The first algorithm introduced in this paper makes extensive 
use of the controllability vector CV. Dealing with vectors 
for the placement procedure allows a better “visualization” 
of the process and, hopefully, allows a more intuitive way 
to explain the obtained results. Further, taking advantage of 
linear algebra properties and metrics, such as orthogonality 
and norms, we are able to save significantly on computational 
effort, thus reducing the computational time. In the following 
sections we outline the algorithm based on definition ([2]). 

A. Placement Criteria 

The placement of the controller should be based on two 
main criteria. First, the infiuence of the controller on a specific 
line. Seeking to maximize controllability as defined in ([S]), the 
most effective controller placement is at a location where a 
marginal change in the controller setpoint will result in the 
maximum change of the line flow. The second criterion is 
the number of individual lines it can control in an effec¬ 
tive manner. In the ideal case, we are seeking a controller 
placement which can control as many lines as possible, and, 
at the same time, it will need the least control actions to 
achieve a power fiow change. As these two objectives cannot 
always be satisfied at the same time, we seek the best possible 
compromise. We will see in the following how we translate 
these two objectives in equivalent metrics, and how we design 
our algorithm in order to identify the optimal controller 
placement. 


B. Metrics for the placement criteria 

Our approach is based on the controllability vector CV in¬ 
troduced in Section |III1 A suitable placement metric would be 
the 1-norm of the controllability vector CV, for two reasons: 
first, it takes into account all vector elements (in comparison 
with the cxo-norm), and second it handles uniformly all vector 
elements, i.e. it does not give an advantage/penalty to extreme 
values, such as the L2 and higher order norms. 

The placement of the second controller should infiuence a 
set of lines as different as possible from the first controller: 
we are thus seeking a CV orthogonal to the first one. Still, the 
vector norm should be taken into account, as this determines 
the degree of infiuence. In this paper, we select as metric the 
volume of the poly tope defined by the vectors in question. 

Here, we will focus on the motivation for using the CV 
vectors and the volume as placement metric. We define the 
poly tope from the projections of the vectors on the coordinate 
axes, as explained in the following. In general, the larger the 
volume covered by the polytope, the higher and more diverse 
is the infiuence of the controllers on the line fiows. This is 
schematically presented in Fig. [U Assume that axes x, y, z 
represent the fiows on three AC lines. The orange and green 
vectors are controllability vectors, representing the infiuence of 
a controller placement on these lines. The higher the projection 
of the vector on each axis is, the higher the effect that a 
marginal change of this controller has on the specific line fiow. 
The convex hull of each CV is determined from the projection 
of the CV on each axis. As we can observe the green vector has 
a substantial infiuence on two of the lines (along x and y axes), 
but it does not affect significantly the third line (along the z 
axis). On the other hand, the orange controllability vector has 
a more balanced infiuence on all lines. Selecting the 1-norm 
as a metric, the green vector has a higher 1-norm than the 
orange vector. On the other hand, the convex hull resulting 
from the orange vector results in a larger volume than the 
convex hull from the green vector and has a more uniform 
effect on a larger set of lines. In Section IVIII-DI we compare 
the two placement metrics. 



Fig. 4. Schematic representation of the 1-norm and the convex hull volume 
metrics. The orange controllability vector has a lower 1-norm but its convex 
hull covers a larger volume than the green vector. 

To calculate the convex hull, we transform the vector 
CV = [ci] into the matrix X = diag{[ci]), where each row 
corresponds to a point in the -dimension space. We augment 
matrix X by adding the point of origin, i.e. a zero vector, so 
that X' = diag{ft Ci]). The set of these points define a conical 
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hull, which is by definition convex. Convexity is a necessary 
condition in order to compute the resulting volume. 

As shown in ([T71) . the conical hull is the set of all conical 
combinations for a given set S. 


coni{S) = 


. i=l 


aiXi\xi e S,ai e M, Qfi > 0,i,/c = 



Based on dU, (|5l), (fTTl) . vectors Xi are the row vectors of 
matrix X', while the scalars represent APdc- Contrary 
to this definition, APdc could also be negative. In order to 
account for that in our algorithm, we consider instead both 
CV and —CV for each placement. A detailed description of 
the algorithm follows in the next sections. 


C. CV-based Placement Algorithm 

In order to have a fair comparison among all controllers, we 
assume the same marginal change in the controller fiow for 
all possible placements: APdc = 1. As placement measure 
we use the sum of these volumes. The larger the sum, the 
more effective the controller placement. For each subsequent 
placement, in order to avoid adding overlapping volumes, we 
determine the extreme values for each orthant and calculate 
the resulting volumes. To reduce the number of subsequent 
placements to be examined, we consider CV which are as 
orthogonal as possible to the vectors we have already selected 
0. As shown in Figure [51 such vectors would not necessarily 
result in the largest possible convex hull volume, but they will 
offer more diverse controllability options. 



Fig. 5. Schematic representation of the convex hull volumes resulting from 
two nearly orthogonal vectors (green and red; their addition is shown in blue) 
in comparison with the volume of the convex hulls resulting from vectors not 
too orthogonal (green and black; their addition is shown in purple). 

The algorithmic steps are outlined below: 

1) Form diagonal matrices from the CV, including the 

point of origin: = [0; dia^([c^’^])] 

2) Calculate V{^^ = 

3) Place the first controller between the nodes m,n that 
resulted in the largest V{^^. 

^By definition, two vectors can either be orthogonal or not. Here, by the 
phrase “as orthogonal as possible” we imply that the angle between the two 
vectors should be as close to 90^ as possible. For the sake of readability in 
this and the following sections, we will abuse the term “orthogonal” in order 
to compare the angles between vectors, e.g. “less orthogonal” will denote that 
the angle between a pair of vectors deviates further away from 90°. 


4) Calculate the orthogonal projections of all other con¬ 
trollability vectors C Vij on the vector space defined by 
CV„,n. 

a) Calculate matrix: Pa = A(A^A)“^A^, where 
A=[CV„,n] 

b) Calculate the orthogonal projections: = 

PaCVij, for i,jf ^ m,n 

5) Calculate: cos(pij = holds cos<pij G [0,1]. 

6) Select the /c-vectors CVij with the lowest cos(j)ij . 

7) For all /c-vectors formulate matrix S = [CVi 

-CVi CV2 -CV2 CV1PCV2 CV1-CV2 

-CVi + CV2 -CVi - CV2]. 

8) Divide the column vectors of matrix S in the orthants 
they belong to. For each orthant p form a matrix 

with the respective column vectors. 

9) Select the extreme values from each row vector of matrix 

gorth,p vector ^ order 

to avoid adding overlapping volumes). 

10) Compute the volume V^^^ of the conical hull defined 

from the row vectors of matrix _ 

nn orth—max,p-\\ 

diag{[0; ’^]). 

11) Repeat Steps [9] and [T0| for all p orthants. 

12) Compute the sum of the conical hull volumes for this 

placement: = J2p ^ 2 .p- 

13) For the second placement, select the nodes j, for which 
V2 becomes maximum. 

14) For subsequent placements go to SteplH now assuming 
the vector space will be defined by the vectors A = 
[CV^n CVij . . .] 

D. Case Study 

a) Comparison of the two metrics: 1-norm volume of 
convex hull: Consider the power system model of Fig. |2l In 
this paragraph we compare the 1-norm metric with the volume 
of the convex hull defined by each CV. We used function 
convhulln of Matlab for the volume calculation. 



Fig. 6. 1-norm of Controllability Vectors CVmn- The vectors are sorted from 
minimum to maximum volume of the convex hull they define. The blue dashed 
line represents the trendline, denoting a correlation between the convex hull 
volumes and || CVmn ||i 

Figure [6] presents the 1-norm of all controllability vectors, 
sorted from minimum to maximum volume of the convex hulls 
they define. The blue dashed line shows the trendline of the 1- 
norm vectors - it is evident that there is a correlation between 
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the two metrics. Pair 4-8 results in the highest “controllability” 
with both metrics. The bottom-five placements are the same for 
both metrics, nevertheless with a slightly different ranking. For 
the reasons outlined above, the volume of convex hull seems to 
represent better the degrees of freedom that are added through 
the placement of each controller. 

b) Placement of three HVDC lines: As shown from 
the results presented in Fig. [6l the first HVDC line should 
be placed between nodes 4 and 8. Figure [7] sorts the 1- 
norm vectors based on the polytope volume for the second 
placement. The CVij “most orthogonal” to CV 48 , i.e. with 
cos(j) < 0.2 are colored red to yellow. The CVu is selected 
for the second placement. We repeat the same procedure for 



Placement of Second HVDC line 


Fig. 7. Second HVDC placement. First HVDC placed between 4-8. The 
CVs are sorted with increasing volume of the convex hull they create in 
combination with CV48- Red/yellow-hued bars denote the most orthogonal 
CVs, i.e cos(j)(^cVij-CV 48 ) — closest to zero, yellow closest to 

0.2. The CV to be selected is one with coscp < 0.2 that results in the largest 
volume. Here, this is between nodes 1-7. 

the third HVDC line placement (figure not shown due to 
space limitations). In this case, only two vectors are below the 
“orthogonality” limit, CVg-io and CV 2 - 10 , thus the algorithm 
saves significantly on computation time. The third line is 
placed between nodes 9 and 10. 

In the following section, we compare these results with 
an optimal placement algorithm we developed based on an 
alternative definition for controllability. 

IX. Optimal Placement based on an Alternative 
Deeinition of Controllability 

A definition for controllability alternative to (O is the 
following: we define controllability as the sum of control 
actions |APc(j)| for changes in line fiows \APL{i)\ < p, on 
all AC lines i G £ = {1,... 

ric 

= ^TOm|APc(i)|, (18) 

i=i 

SO that Vi G £, |AP^(^)| can take any value between 
Assuming the DC approximation for the power flow, 
will be linear in p. Compared with here we define 
controllability as the magnitude of the control action necessary 
to perform a change in each of the line fiows, instead of the 
magnitude of the change in line flows for a marginal controller 
change; so here we are looking for the minimum Based on 


([18]) and vector CV, we formulate the following optimization 
problem. For every different HVDC placement and each “set” 
of AC line fiows: 

minl^t (19) 

subject to: 

-t < < t 

0 < t 

By the term “set”, we denote the set of AC lines to be 
independently controlled during each HVDC placement. For 
example, two HVDC lines can independently control two 
AC lines, so for the second HVDC placement we run the 
algorithm iteratively for all pairs of AC lines and we add the 
magnitude of control actions. We test three different ways to 
determine APl. First, constant APl^i) = lOOMIV for all 
i e riL. Second, APl^j) proportional to the line limit of line 
i (APL(i) = 0.1Pi;,(^)). Third, APl^j) proportional to the line 
reactance of line i (APl^j) = 1000 • with XLi in p.u. 


CV • = APl 

pDC / pDC / pDC 
max — max 


( 20 ) 


100 

80 

60 

40 

) 

20 


-■-ACLim 
- Const 
-•-ACLineX 


mri 

; f f:*f ^ 

' A A A A <b<o<b<b^<o Wb9: 



Placement of Second HVDC line 


Fig. 8. Second HVDC line - alternative controllability definition. Sum of 
control actions for each placement (expressed in % relative to the worst case). 
[ACLim]: APL(i) proportional to the line limit; [Const]: APl the same for 
all lines; [ACLineX]: APl proportional to line reactance. 
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Placement of Third HVDC line 


Fig. 9. Third HVDC line - alternative controllability definition. Sum of 
control actions for each placement (expressed in % relative to the worst case). 
[ACLim]: APL(i) proportional to the line limit; [Const]: APl the same for 
all lines; [ACLineX]: APl proportional to line reactance. 

For the first HVDC placement, the algorithm minimizes 
||ft)c||i- This is heavily infiuenced by ||CV||i which is 


























































presented in Fig. [6l Same to the algorithm presented in 
Section IVIII-Cl this algorithm also results in the selection of 
CV^s as the best placement for the first HVDC line. 

Figure [8] presents the results for the second HVDC place¬ 
ment based on definition ([18]). We observe here that all three 
approaches for setting APl result in similar ranking for the 
top five possible locations. Again, both algorithms result in 
the same placement for the second HVDC line CV 17 . Fig. [9| 
presents the results for the third HVDC placement. As we 
see here, the three different APl approaches result in similar 
ranking for the first nine possible locations. On the other 
hand, we observe that the two placement algorithms, based 
on alternative definitions, result in different optimal placement 
locations for the third line. 

X. Discussion 

Comparing the two placement algorithms, although their 
objective functions are different, they can arrive to similar 
results if the optimal placement fulfills both definitions for 
controllability. The advantages of the vector-based algorithm 
is that it requires substantially less computation time and that 
we do not need to determine different metrics for setting APl. 
The LP-problem, in comparison, requires 4’095 iterations for 
the second placement and 16’380 iterations for the third 
placement. On the other hand, the vector-based algorithm 
requires a proper calibration of the heuristics used (e.g., 
limit). To eliminate redispatching costs with a limited number 
of controllers, besides controllability, loading patterns and 
locations and costs of the generators need to be considered. 
This paper focused on the introduction of the concept and 
the placement algorithms. Future work will extend these 
algorithms, and will also investigate how different APjjc 
setpoints for different lines affect the algorithm performance. 

Placing HVDC lines in parallel: In Fig. [8] and O we 
observe that the magnitude of control actions is zero when 
placing an additional HVDC line in parallel with an existing 
one. Such a placement results to an infeasible solution. This 
is because an arbitrary number of parallel HVDC lines cannot 
control the flow of more than one AC line in the system. With 
APl = CV^n(P^„^(i) +P^f( 2 )) and being 

scalars, by placing two HVDC lines in parallel we maintain 
only one degree of freedom. 

Application for TCSCs and PSTs: For TCSCs, instead 
of the controllability vectors the relationship in (O can be 
used, where and APl = • Axij. 

The difference is that in the case of TCSCs, (fl^ is dependent 
on the operating point - the vector Pb is part of the equation. 
As a result, the placement algorithm might need to take 
several instances of (Hsii into account for each pair of nodes, 
considering different operating points. For PSTs, an expression 
similar to (O can be derived, where the derivative to be 
computed will be with respect to the angle change instead 
of the reactance change. Again, it is expected that this vector 
will depend on the operating point as with TCSCs. 

XL Conclusions and Outlook 

In this paper we introduced the concept “A Fully Con¬ 
trollable Power System”, which aims to fully decouple the 


market operations from the security considerations in power 
systems. In a fully controllable system, the electricity market 
could freely determine the generator dispatch, keeping the 
redispatching actions which maintain system security to a 
minimum. In comparison with current AC power systems, 
controllable line fiows eliminate the dependency on the in¬ 
jection points. Thus, in case of contingency the power can 
be rerouted so that all load is served and no overloadings 
occur. This minimizes the redispatching costs, and in effect 
eliminates the Cost of Security. 

We have dealt with three main topics in this paper. First, 
we determined the upper bound on the number of controllable 
lines in any system, which is ul — riB P I for placing 
controllers in series and — 1 for placing controllers in 
parallel with existing AC lines (jil is the number of lines 
and riB is the number of nodes). Showing that each controller 
can control a single AC line, the same expressions serve as 
lower bounds for the necessary number of controllers to make 
a power system “fully controllable”. Second, we intoduced the 
controllability vector, as a measure of the controller effect on 
any AC line fiow. Third, we used this vector to formulate two 
controller placement algorithms that maximize controllability 
and compare their results. The formulations depend on two 
alternative definitions of controllability. 

Future work includes the formulation and demonstration of 
a convex optimization algorithm based on the controllability 
vector, alternative to the algorithm presented in IVIII-CI 
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